Measuring the transmission phase of a quantum dot in a closed interferometer 
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The electron transmission through a closed Aharonov-Bohm mesoscopic solid-state interferome- 
ter, with a quantum dot (QD) on one of the paths, is calculated exactly for a simple model. Although 
the conductance is an even function of the magnetic flux (due to Onsager's relations), in many cases 
one can use the measured conductance to extract both the amplitude and the phase of the "intrinsic" 
transmission amplitude to = —i\to\c laD through the "bare" QD. We also propose to compare this 
indirect measurement with the (hitherto untested) direct relation sin 2 (ao) = |£d| 2 / max(|£_D| 2 ). 

PACS numbers: 73.63.-b, 03.75.-b, 85.35. Ds 



Recent advances in nanoscience raised interest in quan- 
tum dots (QDs), which represent artificial atoms with 
experimentally controllable properties [1,2]. Connecting 
the QD via metallic leads to electron reservoirs yields 
resonant transmission through the QD, with peaks when- 
ever the Fermi energy in the leads crosses a resonance on 
the QD. The energies of the latter are varied by con- 
trolling the plunger gate voltage on the QD, V. The 
quantum information on the tunneling of an electron is 
contained in the complex transmission amplitude, in = 



The phase od is particularly interesting, 



given its relation to the additional electron occupation in 
the system via the Friedel sum rule [3,4]. This phase is 
also predicted to exhibit interesting behavior e.g. near a 
Kondo-like resonance [5]. This motivated experimental 
attempts to measure an [6,7], using the Aharonov-Bohm 
interferometer (ABI) [8]. 

In the ABI, the QD is placed on one branch, in par- 
allel to a "reference" branch (both connecting the two 
external leads). A magnetic flux $ in the area between 
the two branches creates a phase difference 4> — e&/frc 
between the wave functions in the two branches [9]. In 
the two-slit limit, the total ABI transmission is 

T=\t\ 2 = \t D e^+t B \ 2 =A + Bcos{cj) + (3), (1) 

with [3 = old — k, where k contains V-independent 
contributions from the reference transmission, ts = 
—i\tB\e l5B , and from the electron "optical" paths on the 
two branches. However, for the "closed" two-terminal ge- 
ometry, unitarity (conservation of current) and time re- 
versal symmetry imply the Onsager relations [10], which 
state that the two-terminal conductance, G = (e 2 /h)T, 
is an even function of (p. Therefore, a naive fit of the 
experimental transmission to Eq. (1) must yield /3 = 
or 7r - with no relation to ud- Indeed, the experimental 
data [6] for T depend only on cos</> [11]. 

Aiming to measure a non-trivial AB phase shift f3 then 
led to experiments with "open" interferometers [7,12], 



which contain additional "leaky" channels, breaking the 
Onsager symmetry. A fit to Eq. (1) then yields a phase 
(3 which increases (with V) gradually from to it through 
each resonance. However, the detailed ^-dependence of 
(3 depends on the strength of the coupling to the addi- 
tional terminals [13]. Although it is possible to optimize 
this strength, and reproduce the two-slit conditions [14], 
this involves large uncertainties. 

In the present paper we present exact results for the to- 
tal transmission of the closed ABI, T. Although T is even 
in (p, contradicting the simple two-slit Eq. (1), it does de- 
pend on both Td and arj. Under appropriate conditions 
(see below), one can thus extract an from the measured 
T, eliminating the need to open the interferometer. This 
possible extraction was not noticed in earlier discussions 
of the closed ABI. Theoretical analyses used the Keldysh 
technique, combined with the wide-band and related ap- 
proximations [15,16], or ignored electron-electron interac- 
tions [17]. These approximations, which sometimes miss 
important features of the results (see below) , are avoided 
in our calculation, which is done in the linear response 
limit, and at temperature T = 0. 

We demonstrate our results for a simple lattice model, 
shown in Fig. 1: for $ — 0, each (unit length) segment 
in the figure represents a real tight-binding hopping ma- 
trix element — J, — II, — Ir, — Jl and —Jr, as in- 
dicated. All the on-site energies are zero, except tu on 
the site "dot" and eo on the site "ref" (which sits on 
the reference path and represents a simple point contact, 
a tunnel junction, etc). The latter two energies can be 
varied experimentally by the plunger (or point contact) 
gate voltages eo = V and eo = Vq [11]. As usual for such 
models, electron-electron interactions are included only 
via an on-site Hubbard interaction U on the QD. The AB 
phase in the triangle, <fi = 4>l +<^ii) is included by attach- 
ing a factor e lc ^ L (e % ^ R ) to the hopping matrix element 
Jl {Jr)- At T = 0, the electron energy = — 2 J cos k 
is equal to the Fermi energy on the leads, ej?, and we 
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calculate the transmission for electrons with spin a. 

We start by reviewing the "intrinsic" transmission 
through the QD, without the reference path (e.g. for 
large |Vo| — | e 1 , or with II = Ir = 0). Adapting the 
results of Ref. [18], one has 

tD = -ijDsin-a D e laD = 2ism\k\J L J R g D (e k )/ J, (2) 

with the QD asymmetry factor jd = 2JlJr/(J£ + J R ) 
and the "intrinsic" Green function on the QD, .gzj(efc) — 
l/[efe - £d - T, D (e k )}. Here, Sn(e fc ) is the self-energy 
on the QD, which contains contributions from the leads, 
^D.ext = -e*l fc l(j| + J R )/J (which exists also for the 
non-interacting case [13]), and from the electron-electron 
interactions on the QD itself, Y>D,int(u) (which vanishes 
when U = 0). As £d = V increases, old grows gradually 
from zero (far below the resonance), through n/2 (at the 
resonance), towards 7r (far above the resonance). 

Interestingly, for this one-dimensional model, normal- 
izing the measured To = |tzj| 2 = 7 D sin 2 (a£>) by its 
(V- independent) maximum j D yields the value of a^. 
Assuming coherence, this (hitherto ignored) method for 
measuring old directly from Td eliminates the need for 
any complicated interferometryl [19] In the remainder of 
this paper we discuss ways of extracting old indirectly, 
from the closed ABI measurements. Comparing results 
from sin 2 (a£)) = Td/j%, from the closed ABI (below) 
and from the open ABI [14] (all with the same QD) 
should serve as consistency checks for this conclusion. 

The same analysis yields the transmission amplitude 
through the reference path (when e.g. Jl = Jr = 0), 

t B = -ilB sin 8 B e iSB = 2i sin \k\I L I R g B {e k ) / J (3) 

with the bare reference site Green function g B — l/[e k — 
e + e l|fc| (/I+4)/J], and 7B = 2I L I R / (I 2 L + I R ). In the 
two-slit situation, Eqs. (2) and (3) suffice to determine 
the overall transmission, as in Eq. (1). However, the sit- 
uation is more complicated for the closed ABI. The main 
result of the present paper concerns the exact transmis- 
sion amplitude through the closed ABI, 



t = A D t D e l,l > + A B t B , 



(4) 



where we find A D = g B {ek-^a)G D {e k ) / g D {e k ) and A B = 
l + G D (e k )Y, ext (e k ). Here, G d (uj) = l/[u) - e D - T,(co)] is 
the fully "dressed" Green function on the QD, with the 
dressed self-energy E = Si„ t + Yj ext . Both terms here 
differ from their counterparts in the "intrinsic" £d, by 
contributions due to the reference path. Equation (4) 
looks like the two-slit formula, t — tpe 1 ^ + t B . However, 
each of the terms is now renormalized: Ad contains all 
the additional processes in which the electron "visits" the 
reference site (Ad = 1 when II = Ir = 0), and A B con- 
tains the corrections to t B due to "visits" on the dot. In 
fact, a physical derivation of Eq. (4) amounts to starting 
from Eq. (2), and adding an infinite power series in II 



and Ir. We now discuss the (^-dependence of T = |t| 2 , 
in connection with the Onsager relations and with the 
possible indirect extraction of old- 

We first note that both parts in S(efc) are even in <fi, 
due to additive contributions (with equal amplitudes) 
from clockwise and counterclockwise motions of the elec- 
tron around the ring (see e.g. Refs. [8,13,17,20]). In 
order that T also depends only on cos <p, as required by 
the Onsager relations, the ratio K = A B t B / (Aota) = 
i[G_D(efe) _1 + T, ext (e k )], with the real coefficient x = 
IhlRl\JhJR{tk - eo)L mu st be real, i.e. 



3[Gr)(efc) 1 + T, ext (e k )] = 3£ ir) 



0. 



(5) 



The same relation follows from the unitarity of the 2x2 
scattering matrix of the ring. This relation already ap- 
peared for the special case of single impurity scattering, 
in connection with the Friedel sum rule [4], and was im- 
plicitly contained in Eq. (2), where SS^.mt = [18]. 
Equation (5) implies that (at T = and oj = e k ) the 
interaction self-energy T,i nt (e k ) is real, and therefore the 
width of the resonance, 3G , £)(efe)~ 1 , is fully determined 
by the non-interacting self- energy 9E ext (efe). 

Since T, ext (uj) depends only on the (non-interacting) 
tight-binding terms, it is easy to calculate it explicitly. 
We find T, ext (e k ) = S 

D,ext 

(e k ) + A ext , where 

A ext - e 2 ^g B ( J 2 L ll + J 2 / 2 + 2J l J r I l Ir cos </>)/ J 2 . (6) 

The term proportional to cos0 comes from the elec- 
tron clock- and counterclockwise motion around the ABI 
"ring". Similarly, one can write S mt (e fe ) = £r>,mt(e/s) + 
Aj„ t , and thus Gd^) -1 = gD^k)^ 1 - A, with A = 
A ext + A int . Hence, t = A D t D (e i ' t ' + K). Writing also 
Ad = C/[l - g D (e k )A], with C = (e k - e )g Bl we find 



T= \C\ A T D - 



1 + K 2 + 2K cos <j) 
- 2^[g D A] + \g D A\< 



(7) 



Equation (7) presents an alternative form of our main 
result. Although the numerator looks like the two-slit 
Eq. (1), with [3 = or n (depending on signif), the 
new physics is contained in the denominator - which be- 
comes important in the vicinity of a resonance. The cen- 
tral term in this denominator depends explicitly on the 
phase of the complex number go- Since this number is 
directly related to to, via Eq. (2), one may expect to 
extract old from a fit to Eq. (7), taking advantage of the 
dependence of the denominator on cos <p. Physically, this 
dependence originates from the infinite sum over electron 
paths which circulate the ABI ring. The rest of this pa- 
per is devoted to the conditions for such an extraction. 
Generally, this is not trivial, as one needs the detailed 
dependence of A on cos <p and on the various parameters. 
We have presented this dependence for A ext , but not for 

Aint- 
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The extraction of etc becomes easy when one may 
neglect A int . The simplest case for this is for single- 
electron scattering, when T, int — 0. Interactions (i.e. U) 
are also negligible for a relatively open dot, with small 
barriers at its contacts with the leads [21]. Another 
effectively single-electron scattering case arises near a 
Coulomb blockade resonance, when the effect of inter- 
actions can simply be absorbed into a Hartree-likc shift, 
e_D + Sj„ t — > cd + NU , if one assumes that N depends 
smoothly on the number of electrons on the QD, and not 
on <f> [20]. If one may neglect Ai nt , then A w A ext is 
given in Eq. (6). Using also Eqs. (2) and (3), we find 



r = \c\ 2 t d 



1 + K 2 + 2Kcos< 



1 + 2P(z + cos <j>) + Q(z + cos (j)) 2 



(8) 



where z = (J 2 L I 2 L + J 2 R I 2 R )/(2J L J R I L I R ), P = X[vt B t D ], 
Q = |^b| 2 Td, and v = e 2% ^/{2 sin 2 |fc|) depends only 
on the Fermi wavevector k, independent of any detail of 
the ABI. A 5-paramctcr fit to the explicit 0-dcpendence 
in Eq. (8) for given values of V and Vo then yields 
\C\ 2 T D , K, z, P and Q, and thus cos(a D + S B + 2|&|) = 
P/^/Q, from which one can extract the U-dependence 
of old ■ The same U-dependence of old is also contained 
in K oc (cot art + cot |fc|)). As discussed after Eq. (2), 
our model also implies that To = 7^ sin 2 (an)- Since the 
U-dependence of To can also be extracted from the fit- 
ted values of either \C\ 2 Td or Q, we end up with several 
consistency checks for the determination of old- Addi- 
tional checks arise from direct measurements of Td and 
T R = \t B \ 2 , by taking the limits \Vq\ = | e 1 — ► 00 or 
\V\ = \e D \^co. 

The LHS frame in Fig. 2 shows an example of the V- 
and (^-dependence of T for this limit (no interactions), 
with k — tt/2 and J L — J R = I L = I R = 1, Vo = 4 (in 
units of J), implying K = eo/eo = V/Vq. Far away from 
the resonance T < 1, Q < |P| < 1 and \K\ > 1, yield- 
ing the two-slit-like form T w A + B cos <j>, dominated 
by its first harmonic, with B/A w 2\K~ X — P]. How- 
ever, close the the resonance T shows a rich structure; 
the denominator in Eq. (8) generates higher harmonics, 
and the two-slit formula is completely wrong. This rich 
structure may be missed if one neglects parts of the (jy- 
dependence of A, as done in Ref. [16]. Note also the Fano 
vanishing [22] of T for V <~ 10 at <f> = 2nir, with integer 
n. Without interactions, everything can be extended to 
a QD with many resonances, e.g. due to Coulomb block- 
ade shifts in the effective td with the number of electrons. 
Using a generalization to Eq. (8), given in Ref. [14], the 
RHS frame in Fig. 2 shows results for two resonances, 
with en — ±5. Interestingly, Fig. 2 is qualitatively simi- 
lar to the experimentally measured transmission in Ref. 
[11]. However, so far there has been no quantitative anal- 
ysis of the experimental data. 

To treat the general case, we need information on A int . 
First of all, we emphasize that a successful fit to Eq. 



(8) justifies the neglect of the ^-dependence of Ai nt . If 
the various procedures to determine a B from Eq. (8) 
yield the same F-dependence, this would also confirm 
that Ai n t is negligibly small. A failure of this check, or 
a more complicated dependence of the measured T on 
cos 4>, would imply that A int is not negligible. 

As seen from Eq. (6), A ext is fully determined by a 
single "visit" of the electron at "ref. For small Tb, or 
large \Vq\ = |eo|, it is reasonable to conjecture that A int 
is also dominated by such processes. In that case, we ex- 
pect Ai nt to be proportional to the same brackets as in 
Eq. (6), i.e. A int w w(z + cos^), with a real coefficient 
w. This yields the same dependence of T on cos <f> as in 
Eq. (8), with a shifted coefficient v. If w depends only 
weakly on V, then this shift has little effect on the de- 
termination of a B ■ Again, the validity of this approach 
relies on getting the same U-dependcnce of an from all 
of its different determinations. 

The situation becomes more complicated near a 
Kondo-likc resonance. Maintaining the (non-trivial) as- 
sumption that Gd = l/[ w — eu — S_d(w)], the Kondo 
peak at the Fermi energy must be generated by Ed- For 
the intrinsic QD, this yields old = tt/2 and to = Jd, 
resulting in a V- independent plateau for To . A priori, 
it is not obvious what happens in the presence of the 
"reference" path. Hofstcttcr et al. identified the Kondo 
region by requiring that the phase <5 res of the fully dressed 
Green function Gd be equal to tt/2. Our result for Gd 
shows that this is impossible: the phase <5 ros depends on 
4>, via the (^-dependence of A, and thus cannot be set at 
the constant value tt/2. (Apparently, this (^-dependence 
was neglected in the analytic parts, and weak for the nu- 
merical parameters used in Ref. [16]). Alternatively, one 
might assume that the "bare" QD sticks to the Kondo 
resonance, and thus qc = tt/2 (independent of V) even in 
the ABI. Equation (7) then replaces the Kondo plateau 
by a complicated dependence on <p (including the first 
harmonic), which differs significantly from that of Ref. 
[16]. Clearly, this limit requires more research. 

Finally, we give some more details of our derivation. 
Our Hamiltonian, which simply adds the reference path 
to that of Ng and Lee [18], is 

H = e D ^2 d\d a + — ^ n-danda 

(7 <J 

ka ka 

+ £ o c o<t c o° + ( u k c L c ka + K c L c ov) . ( 9 ) 

o kcr 

where c\ a creates single particle eigenstates (with spin 
cr) on the unperturbed "background" chain (with II = 
I R = J, Jl = Jr = 0), with cigenenergy efe = — 2Jcosfc, 
while c a = Ek cka/VN, U k = -[(/£ - J)e~ ik + [I R - 



J)e ik ]/VN, and V k = -{J L e 



i(j>p—ik 



J Re' 



-i(f> r -\-ik 



)/VN. 
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The operators on the dot, d a and d^, anti-commute with 
Cka, 4o- Also, n da = d\ d a , and ct = -a. 

As stated above, one can derive Eq. (4) by expand- 
ing Eq. (2) in powers of II and Ir. A more general 
approach uses the standard relation between the 2x2 
scattering matrix Tkw and the matrix of retarded single- 
particle Green functions, GWM = <Ws° + 9k T kk'9k'^ 
with <?5?(w) = l/(w — efe), evaluated on the energy shell, 
u = cf = e/s = £fe' [5]. The equation-of-motion (EOM) 
method is then used to express (ui — ek)Gkk>(u) and 
{uo — ek)Gkd{u) as linear combinations of each other and 
of G_d(oj), allowing us to express each of them (and thus 
also t cx Ti fe i i fe i) in terms of G_d(oj), yielding Eq. (4). 
Since we do not use an explicit solution for Gd{u) itself, 
we don't need to deal with the higher order correlation 
functions (due to U), which appear in its EOM. 

We hope that our paper will stimulate attempts to fit 
experimental data to our Eq. (8), and to compare the re- 
sulting o>d with its direct estimate via T. This procedure 
should work in many cases. We also hope that our pa- 
per will stimulate more detailed theoretical calculations 
of Aj n (. As explained, the existing approximate calcula- 
tions miss the crucial 0-dependence of these interaction- 
dependent terms. 
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FIG. 1. Model for the closed ABI. 




FIG. 2. AB transmission T versus the AB phase <f> and the 
gate voltage V, for one (LHS) and two (RHS) non-interacting 
resonances. 
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